Exact dark state solutions of the coupled atomic-molecular Bose-Einstein condensates 

in an external potential 



O 

(N 

as: 
< 

in 



C/2 

c3 



i 

c 

o 
o 



> 
o 
o 
o 
in 

00 

o 



X 



Xiao-Fei Zhang 1 ' 3 , Jun-Chao Chen 2 , Biao Li 2 , Lin Wen 1 , and W. M. Liu 1 

1 Beijing National Laboratory for Condensed Matter Physics, 
Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 
2 Nonlinear Science Center, Ningbo University, Ningbo 315211, China and 
3 College of Science, Honghe University, Mengzi 661100, China 
(Dated: August 26, 2011) 

We consider a coupled nonlinear Schrodinger equations describing an atomic Bose-Einstein con- 
densates coupled to a molecular condensates through the stimulated Raman adiabatic passage loaded 
in an external potential. The existence of dark state are investigated within the full parameter space 
accounts for all the nonlinear collisions, together with the atom-molecule conversion coupling and 
external potential. The results show that there are a class of external potentials such as double-well, 
periodical and double periodical, and harmonic potentials, where the exact dark solutions can be 
formed. Our results may raise the possibility of relative experiments for dark state in the coupled 
atomic-molecular Bose-Einstein condensates. 
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I. INTRODUCTION 

The experimental realization of Bose-Einstein conden- 
sates (BECs) (H-H arouse great interest in ultracold 
molecule due to their potential applications for the tests 
of fundamental physics and for the drifts of fundamen- 
tal constants [j-Q ■ This fundamental system is of spe- 
cial interest since the collisional coherence can led to re- 
markable demonstration of the reversible atom-molecule 
formation thus open the floodgate for many beau- 
tiful experiments, such as coherent oscillations between 
atoms and molecules, dissociation of molecular BECs, the 
"clumping" of the condensates due to the modulational 
instability, and so on [Tol - tl^ . It also represents a rich 
ground for investigations of the nonlinear excitation of 
the coupled atomic-molecular Bose-Einstein condensates 
(AMBEC) [3. 

In real experiments, the cold molecules can be pro- 
duced from a Fermi gas of atoms EMI or an atomic 
BEC based on Feshbach resonances, Raman photoasso- 
ciation, or stimulated Raman adiabatic passage (STI- 
RAP) |17H25| . However, among them, the Feshbach cou- 
pling mechanism is restricted to the creation of molecules 
in the highest rovibrational level and is only practica- 
ble for a limited number of systems |26l - |28| . and the 
losses caused by inelastic atom-molecule collisions oc- 
curs at a significant rate. In the Raman photoassociation 
process, the effective conversion rate of a pair of atoms 
to ground molecule is limited by spontaneous emission 
from the excited molecular state; while the STIRAP is 
known to have the highest rate of efficiency when convert- 
ing an atomic condensates into a molecular one based 
on a Raman transition, where the input Raman laser 
pulse couples the molecular levels and reduces sponta- 
neous emission. During this process, the molecular state 
is in a dark superposition state (or coherent population 
trapping state, CPT state), which decouples from the 
light and thus suppresses losses due to spontaneous light 



scattering [29(. In this situation, the excited molecular 
state of this coupled system remains almost unpopulated. 
Thus it is not difficult to see the existence of dark is the 
key ingredient in effective production of molecule through 
STIRAP technique. 



Generally speaking, the dark state cannot exist 
with any external potential as the coupled nonlinear 
Schrodinger equations (see Eq. ^ below) is noninte- 
grable [30h35| . The main purpose of the present paper 
is to find the exact dark state solutions of this coupled 
system with tunable nonlinear interactions and external 
potential. Based on the framework of a three-component 
mean-field model that takes into account all types of 
the mean-field nonlinear atomic and molecular interac- 
tions and external potential, we present a comprehen- 
sive analysis of the dark states in the parametrically cou- 
pled atomic-molecular Bose-Einstein condensates. The 
results show that there are a class of external potentials, 
such as double- well, periodical and double periodical, and 
harmonic potentials, which can be used to support the 
existence of the dark states. 



This paper is organized as follows. In sec. II, we briefly 
outline the theoretical model under study and the un- 
derlying mean-field coupled equations. Then a class of 
simplification is presented and the reduced coupled equa- 
tions, focusing on the CPT solutions, are obtained. In 
Sec. Ill, the analytical methods for solving the coupled 
Schrodinger equations are introduced, and we present 
several classes of exact dark state solutions where dif- 
ferent potentials can be applied to trap the coupled con- 
densates and relationship between system parameters are 
established. A simple summary is given in the final sec- 
tion. 
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II. THE THEORY MODEL 

To begin with, we consider a quasi-one-dimensional ge- 
ometry of the condensates where the transverse traps are 
tight enough, thus the transverse motion of the conden- 
sates are frozen to the ground state of the transverse 
harmonic trapping potential. The coupled Hciscnberg 
equations of motion for the atomic and molecular com- 
ponents are govern by the system of equations as follows: 
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here we define the order parameter ij) a represents an 
atomic species of mass m in a potential Vx(x); while 
?/>2 and -03 represents the ground and excited molecu- 
lar species with a potential V 2 (x) and V 3 (x), respec- 
tively. The nonlinear atom-atom interaction Uxx = 
4irh 2 axx/m, the atom-molecule interaction Ux 2 = U 2 x = 
(3ai2/4an)C/n, and the ground molecule-molecule inter- 
action U 22 = (a 22 /2axx)Uxx, with ay being the corre- 
sponding s-wave scattering lengths modulated by a Fes- 
hbach resonance. Sli is the Rabi frequency correspond- 
ing to the transitions between the atomic state and the 
excited molecular state, and VL 2 is the one between the 
molecular ground and excited states. Finally, in this 
paper, we will only focus on the degenerate case with 
V 2 = V 3 = 2Vx and ignore loss and growth mechanisms 
[35l HH; thus the total number of atomic particles, in- 
cluding pairs of atoms inside the diatomic molecules is 
conserved in this model. 

From a general mathematical and nonlinear physical 
point of view, Eq. ([T]) is an example of a classically non- 
integrable field theory, which needs to be treated numer- 
ically. To obtain the exact dark state solutions, we make 
a class of simplifications. To this end, we first neglect 
the excited molecular state as it has a small population 
compared with atomic and ground molecular states, and 
then set U 3 i = 0. Similar to homogeneous condensates 
cases, the STIRAP model considered in this paper is also 
found to support a CPT state (with ip 3 (x,t) — 0) in the 
form of 



Mx,t) =0i(x) e - v *,V2(M) =0 2 (x)e- 12 ^. (2) 

Normalizing the time and length in Eq. ([1]) by t' = 
w 2 t/2 and x' = ^Jujmjhx (the tilde is omitted for sim- 
plicity in the following discussions). Then inserting this 
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FIG. 1: (Color online) The stationary dark state wave func- 
tions of the atomic and molecular condensates 4>\\(x) and 
021 (x) given by Eq. (|10[) . and the profile of the external poten- 
tial in the limit k — and ao — 0. The parameters are given 
as follows: S — l,ao = 0, Si = ai = Ao = fi = A = /J, = 1, 
and (722 = 1 in Fig. la and (722 = 34 in Fig. lb. As shown in 
this figure, different external potentials but with similar pe- 
riodical structure can support the same stationary dark state 
wave functions of the atomic and molecular condensates. 



equation into Eq. ([T]), one can readily derive the follow- 
ing stationary equations: 



M0i = + (ffn^i + 512^2)^1 + v { x )<t>i> 

1 j2j>. ^ 

t(W2+3210i)02+^Oz#2, (3) 



M02 



4 dx 2 



with 4> 2 = —fl<p 2 /2. Here we use the fact n 2 = 2/ii = 2/x 
and the definition O — fix/ ^2- In Ref. [35|], Eq. ^ is 
reduced to a single ordinary differential equation and the 
possibilities of existence of inhomogeneous dark states 
of atomic-molecular Bose-Einstein condensates loaded in 
special external potential is discussed. In what follows 
we find that the exact dark state solutions can also exist 
in other types of external potentials, such as double- well, 
periodical and double periodical potentials. 



III. EXACT DARK STATE SOLUTIONS AND 
DISCUSSION 

In this section, we introduce a analytical method to 
derive a class of exact dark state solutions for Eq. 
Without a trap, i.e., Vi = 0, the system is homo gene ous 
and the dark state is well studied (see Refs. [30I l3l[ |37| 
and references therein). In the presence of external po- 
tential, we find that some special types of dark state so- 
lutions can be derived via suitable selections of the non- 
linear parameters and external potential. 

To get the analytical dark state solutions, we assume 
the stationary dark state solutions have the following 
forms: 
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r B + Bicn 2 (Aa;,fc)' v ' 

where en is the Jacobi elliptic function and < k < 1 
is the module of Jacobi elliptic functions, Bo ^ —B\. 
At this point, the analytical dark state solutions for the 
atomic and molecular condensates can be written as 



FIG. 2: (Color online) The stationary dark state wave func- 
tions of the atomic and molecular condensates <j>n(x) and 
02i {x) given by Eq. (|10[) . and the profile of external potential 
in the limit k = and ao = Bq = 0. The other parameters 
are the same as in Fig (ft} except #22 = 8. In this case, the 
shape of external potential shows double-periodical structure. 
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here we set ao = 0. In this situation, the external poten- 
tial can be obtained by solving Eq. ^ and reads 



4>i{x) = S\/a + aip, 
4>2{x) = --Cl(ao + ai<p), 



(4) 



with auxiliary function p — p(x) satisfying the following 
equations 



{ dx' 
dx 2 



= ho + hup + h 2 p 2 + h 3 p 3 + h 4 p , 

h\ 3 9 n 3 

= -y + h 2 p + -h 3 p + 2h 4 p , 



(5) 



where 5 = ±1 and ao, 01, ho, hi, h%, /13, h 4 are constants 
to be determined. Now we shift our attention to the ex- 
ternal potential and further assume the external potential 
is in the form of 



V = 



b + bi<p + b 2 p 2 + b 3 p 3 + b 4 p 4 
d + dip 



(6) 



with 60, 6 l5 62, 6 3 , 64, do, di are constants to be deter- 
mined. Substituting Eqs. flU) and fl§J) into Eq. ([3]) and us- 
ing the auxiliary Eq. ((SJ) , we can obtain a set of ordinary 
differential equations (ODEs) with respect to the con- 
trollable system parameters {fi, gn, gi 2 , g 22 }. As shown 
in Eqs. (j4|) and ([5]), depending on different choices of the 
form of ip(x), we can get various types of external poten- 
tial and exact dark state solutions through Eqs. ((4]) and 
. It should be mentioned that in order to simplify the 
calculation, we further set S = 1 in the following discus- 
sions. Finally, solving these ODEs, we can obtain two 
families of analytical solutions of Eq. ^ . 



A. The Jacobi cn-wave solutions 

We begin with the Jacobi cn- function solution of p{x). 
In this case, we can get 



Vi(x) = {]^ 2 gi 2 ~n 2 g 22 )pl{x) + {Zgi r 



■ 2512V1 (x) 



h 2 



(9) 

It is evident that the external potential is a function of 
nonlinear parameters, chemical potential, the Rabi fre- 
quencies, and the auxiliary function p{x). In the fol- 
lowing discussions, we find that with different choices of 
the auxiliary function p(x), there exist a class of external 
potentials, such as double- well, periodical and double pe- 
riodical, and harmonic potentials, where analytical dark 
state can exist and the constraint conditions for the non- 
linear parameters can be obtained. For simplicity, we 
only discuss the exact dark state solutions of Eq. © 
in two cases with the module of Jacobi elliptic functions 
k = and k = 1. 

Case 1. In the limit k = 0, one can easily get 
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y B + Bi cos 2 (Ax)' 

— flaiAo 
2B + 2Bi cos 2 (Ax)' 



(10) 



where A = ±1. In this case, the relationship between 
interaction parameters read 
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with Ao = 1 and 522 a function of the system parameters. 

Shown in Fig. [1] are the profiles of the wave func- 
tions for the atomic and molecular condensates </>i {x) and 
4> 2 {x), and the external potential in the limit k = 0. The 
parameters are given as 8 = 1, Bi = ai = Ao = f2 = 1, 
A = 1, and g 22 = 1 in Fig. 1(a) and g 22 = 34 in Fig. 
1(b). According to Eq. (fTTj) . which gives the relation- 
ship among these three kinds of nonlinear interactions; 
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FIG. 3: (Color online) The stationary dark state wave func- 
tions of the atomic and molecular condensates 4>i2{x) and 
4>22{x) given by Eq. (|12[l . and the profile of external poten- 
tial in the limit k = 1. The parameters are given as follows: 
<5 = 1, B = Bi = a = 1, A = 1, ft = -1, and 322 = -20/3 
in (a) and 522 = — 18 in (b). 



for Fig. 1(a), both the atom-atom and atom- molecule 
interactions gu — —2.25 < and 512 = — 7.5 < 0, cor- 
responding to attractive atom-atom and atom-molecule 
interactions; while for Fig. 1(b), gu > and 512 > 0, 
which together with 322 = 34, corresponding to repul- 
sive atom-atom, atom-molecular, and molecule-molecule 
interactions. 

As shown in this figure, different external potentials 
but with similar periodical structure can support the 
same stationary dark state solution. However, we want to 
point out the difference between Fig. [IJa) and Fig. [TJb). 
For Fig. [Ua), * ne external potential is relative shallow 
and the density distribution is mainly locating at the 
crest, which is the result of the competition between the 
nonlinear interaction and the external potential; while 
for Fig. QJb), when the atom-atom and atom- molecule 
interaction become repulsive, the height of the external 
potential increases, and the density distributions of the 
order parameters arc different from the former one and 
mainly locating at the wave trough. 

Figure [5] shows the profiles of the atomic and molec- 
ular condensates 4>\(x) and 4>2{ x ) an d external poten- 
tial in the case of 322 =8. In this case, gu = —0.5 
and <7i2 = —4, which are smaller than the parameters 
used in Fig. 1(a), but with a stronger molecule-molecule 
interactions. It is easy to see that the external po- 
tential presented here shows double-periodical structure, 
which can be realized by two pairs of counter-propagation 
lasers, and widely used in today's ultracold atom exper- 
iments. Finally, no matter the periodical or the double- 
periodical potentials, the density distribution always con- 
tinuous distribution within the extension of the external 
potential. This is different from the localized solutions 
reported in [35| and the following solutions, where the 




FIG. 4: (Color online) The stationary dark state wave func- 
tions of the atomic and molecular condensates (j>i3{ x ) and 
(j>2z{x) given by Eq. (|16[) . and the profile of the external po- 
tential given by Eq. (|17|) . The parameters are given as follows: 
5 = 1, A = 1, ao = 0, ft = 1, A = l,gri2 = l,ai = 1,322 = 
20,n = -1. 



chemical potential must be smaller than zero. 

Thus we conclude that when the ratio of the Rabi fre- 
quency is fixed, both periodical and double-periodical po- 
tentials can be applied to the coupled atomic-molecular 
Bose-Einstein condensates to get the analytical dark 
state solution via varying the strength of the nonlinear 
interactions. 

Case 2. Now we shift our attention to the case of fc = 1. 
In this situation, the stationary solutions for the atomic 
and molecular condensates can be written as 



/ a £?iscch (Aa;) 
V B + Bisech 2 (Aa;) 

— f2ao-Bisech 2 (A:r) 
2B + 2Sisech 2 (A.x) 



(12) 



The relationships among the interaction parameters in 
this case read 
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Shown in Fig. [3] are the wave functions of the atomic 
and molecular condensates <pi2(x) and 022 {%), and the 
profile of external potential in the case of attractive 
molecule-molecule interaction. The parameters are given 
as follows: 5 = 1, Bq = B\ = ao = 1, X = 1, Q = — 1. 
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In the case of 5 = \ = SI = Aq = g\ 2 = 1 , clq = 0,ai 
1 . g 2 2 = 20, /i = — 1, the wave functions of the atomic 
and molecular condensates can be expressed as 



csgn(l / cosh(x)) 

>13 = — ,023 

cosh(x) 



2 cosh(x) 2 ' 



(16) 



with csgn the sign function for real and complex expres- 
sions. In this situation, the external potential is in the 
following form 



FIG. 5: (Color online) The same as in Fig. [4] but f° r a i = 3. 

The molecule-molecule interaction 1722 = —20/3 for Fig. 
in^a) and (722 = — 18 for Fig. Hfb), corresponding to the 
two different cases (i) and (ii) in Eq. (fT4]) , respectively. 
According to the constrain conditions for the nonlinear 
parameters, we can get gu = 5/6, 512 = 14/3, and 
1722 = —20/3 for Fig. EJa), which corresponding to re- 
pulsive atom-atom and atom-molecular interactions and 
attractive molecule-molecule interaction; while for Fig. 
[31(b), the nonlinear parameters read gn = —2, 912 = —1, 
and 522 = —18, where all the interactions are attractive. 
As shown in Fig. [U[a), the external potential shows the 
harmonic structure, while double-well structure in Fig. 
|3Kb), which is similar to the situation in Fig. (fT]): the 
order parameters $12 (x) and (^22 (x) are with the same 
distributions in different potentials. 

It is necessary to point out the stationary dark state 
solutions obtained in this case is different from the ones 
in Case 1 as </>i2(x) and </>22 (x) are localized and sym- 
metric with respect to certain point X (here we choose 
X=0). Thus we conclude that the exact localized dark 
state solutions can be achieved with a proper choice of 
the nonlinear interaction parameters and external poten- 
tial, and both harmonic and double- well potentials can be 
used to trap the coupled atomic-molecular Bose-Einstein 
condensates to get the analytical dark state solution via 
varying the strength of the nonlinear interactions. 

B. The Jacobi sn- and dn-wave solutions 

In this subsection, we consider the Jacobi sn- and dn- 
function solution of <p(x). In this case, one can easily get 



V 2 (x) = 



-1 - 8cosh(x) 5 
4cosh(x) 4 



(17) 



¥2 = A + 



Aisn(Ax, fc) 5 
dn(Ax, k) 2 + 1 



(15) 



where sn and dn arc the Jacobi elliptic functions and 
< k < 1 is the module of Jacobi elliptic functions. 
With the same procedure in subsection A, we can obtain 
a family of exact dark state solutions. As the expressions 
of wavefunctions for atomic and molecular condensates 
are too complicated, here we only focus on the limit k = 1 
and list some main results. 



Figure|4]exhibits the profiles of the order parameters of 
the atomic and molecular condensates </>i3 (x) and (f>23 (x) 
and the external potential. As shown in this figure, the 
external potential is similar to the usual harmonic one 
and the result is similar to the one shown in Fig. (lc) 
in [35j ; while Fig. [5] shows the same wave functions with 
the same parameters used in Fig. [4] but for a% = 3. It 
is easy to see that the external potential in Fig. [5] is dif- 
ferent form the usual harmonic one, but with a peak in 
each side of the potential. It should be mentioned that in 
this case the nonlinear parameters arc different. Thus we 
conclude that by manipulation the nonlinear parameters, 
there exist a class of external potentials where analytical 
dark solution can be obtained. Finally, we want to point 
out that the localized stationary dark state solutions and 
the corresponding two types of external potentials only 
exist for \i < 0, which agrees well with the results ob- 
tained in [35| . 

The stability of the localized stationary solutions ob- 
tained in this paper is still an open problem. This can be 
achieved by carrying out numerical simulations of the lo- 
calized stationary solutions with the imposed initial per- 
turbations, or equivalently adopting the linear stability 
analysis: a small deviation away from the intended state 
is added to see the dynamical evolution of the system; if 
the deviation remains small, then the stability (or adia- 
baticity) is obeyed. Since the object of this work is to 
obtain some exact analytical solutions, a thorough anal- 
ysis of the stability is beyond the scope of the present 
paper and these works belong to further studies. 



IV. CONCLUSIONS 

In summary, we focus on an atomic Bosc-Einstein con- 
densates coupled to a molecular condensates through the 
stimulated Raman adiabatic passage loaded in an ex- 
ternal potential, with the emphasis on the stationary 
dark state solutions. Within the full parameter space 
accounts for atom-atom, atom-molecule and molecule- 
molecule collisions, together with the atom-molecule con- 
version coupling and the external potential, we detailed 
investigate under which conditions the exact dark state 
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solution can exist. Our results show that via suitable 
selections of the nonlinear parameter and external po- 
tential, the analytical dark state solutions can be con- 
structed within a class of external potentials including 
double-well, periodical and double periodical, and har- 
monic potentials. The obtained results are of particular 
significance to experimental realization of such dark state 
in coupled atomic-molecular Bose-Einstein condensates. 
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